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Abstract 

In many cases the correct theoretical description of flexoelectricity 
requires the consideration of the finite size of a body and is reduced to 
the solution of boundary problems for partial differential equations. 
Generally speaking, in this case one should solve jointly the equations 
of polarization equilibrium and equations of elastic equilibrium. How¬ 
ever, due to the fact that typically flexoelectric moduli are very small, 
usually one can consider the solution of polarization equilibrium equa¬ 
tions at a given elastic strain (direct flexoelectric effect) or the solution 
of elastic equilibrium equations at given polarization (converse flexo¬ 
electric effect). Derivation of the polarization equilibrium equations 
and boundary conditions for them can be made in the quite usual way. 
Solution of these equations usually is not too difficult problem. On 
the contrary description of converse flexoelectric effect is more compli¬ 
cated problem. Inter alia effective solution of corresponded boundary 
problems requires the development of special mathematical methods. 
The subject of this paper is a detailed discussion of the relevant theory 
focuses on the converse flexoelectric effect. It is also considered some 
particular examples illustrating the application of the general theory. 


1 Introduction 

A body in which flexoelectrical phenomena are observed always has a finite 
size. In some cases the finite sample size is not significant, and the sample 
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can be described as an infinite medium. Examples are the influence of the 
flexoelectric effect on the dispersion of acoustic phonons [I] or its influence on 
structure of ferroelectric domain wall [2]. In other cases, such as flexoelectric 
bending of a plate [3j, finite size of a sample plays a fundamental role. 

Obviously, the effects associated with the finite size of the sample caused 
by the presence of the sample boundary. In particular, this raises the surface 
piezoelectricity which, as it turns out, leads to the effects of the same order 
as the flexoelectric effect [3]. 

Less trivial fact is that even without surface piezoelectricity, in the pres¬ 
ence of flexoelectricity classical theory of elasticity should be modified. In this 
case it turns out that, generally speaking, one should use elastic boundary 
conditions of non-classical form. These boundary conditions are consistent 
only if one considers the spatial dispersion of the elastic moduli (contribution 
to the thermodynamic potential which is a bilinear in elastic strain gradients) 
[5]. The latter leads to that not only the boundary conditions, but also dif¬ 
ferential equations of elastic equilibrium should be non-classical. 

Equations of elastic equilibrium are needed to describe the converse flex¬ 
oelectric effect (mechanical response to the polarization). As regards the 
direct flexoelectric effect (polarization response to an inhomogeneous elastic 
strain), to describe it the equations of polarization equilibrium are needed. It 
is essential that this equations and boundary conditions for them are modified 
by flexoelectricity not so radically as the equations and boundary conditions 
describing the elastic response (converse flexoelectric effect). This is related 
to the fact that the flexoelectric contributions to the thermodynamic poten¬ 
tial contains only the first spatial derivatives of the polarization, while it 
contains the second spatial derivatives of the elastic displacements. Equa¬ 
tions of polarization equilibrium and boundary conditions for them in the 
presence of flexoelectricity were derived in paper [6], 

Thus, non-trivial features of flexoelectricity in finite samples reveal them¬ 
selves mainly in the converse flexoelectric effect. This is why we do not discuss 
much equations of polarization equilibrium, focusing mainly on the converse 
flexoelectric effect at a given polarization. Moreover, in the particular ex¬ 
amples (but not in the general theory) further we assume a homogeneous 
polarization at which non-trivial features of flexoelectric effect in finite sam¬ 
ples reveal itself most dramatically. 

The case of homogeneous polarization is interesting in the sense that 
a naive analysis based on constitutive equations leads to the (erroneous) 
conclusion that the homogeneous polarization does not cause a flexoelectric 
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deformation of a body. The consequence of this conclusion is the ability 
to create a sensor which not behave as an actuator um- Such an ability 
contradicts the general principles of thermodynamics. But in reality viola¬ 
tions of the principles of thermodynamics is not happening, homogeneous 
polarization does lead to a flexoelectric deformation of the body @1 EJ- 

In the paper [3] the conclusion that homogeneous polarization leads to a 
flexoelectric deformation of a thin plate was made on the basis of the partic¬ 
ular ansatz for elastic displacement distribution over the sample and direct 
minimization of the thermodynamic potential. A more rigorous analysis of 
this problem requires the solution of differential equations of clastic equilib¬ 
rium with appropriate boundary conditions. Such a solution for the case of 
homogeneously polarized ball was presented in the paper Pj. As was ex¬ 
pected the solution confirmed the presence of flexoelectric strains caused by 
homogeneous polarization. 

It is worth mentioning that in the presence of flexoelectricity the exact 
(within the framework of a continuum media theory) equations of elastic 
equilibrium are too complex to be used for a case differs from a simple case 
of a ball. Even for a ball the solution is very cumbersome and requires an 
introduction of non-standard functions presented by a series in powers of 
radial coordinate. So that the development of approximate method is worth. 
Such a method was developed in the paper [10]. 

Papers on the subject 0 [ID] are too brief. More detailed description 
is presented here. General aspects of derivation of the equations of elas¬ 
tic equilibrium and the boundary conditions for them in the case when the 
thermodynamic potential depends on higher derivatives discussed in another 
paper [12]. Here we apply these results to the case of flexoelectricity and 
develop the theory further. Besides, it is presented a comparison of exact 
solution [19] and the solution of the same problem within the approximate 
theory [10] . Some additional examples of approximate theory application 
are presented also. 
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2 Thermodynamic potential, differential 
equations of elastic equilibrium, and 
boundary condition 

It is well known m that within the theory of continuum flexoelectricity is 
described by a contribution to thermodynamic potential density: 

'Hfix = -fj&Pkikjji ~ fSij P k,iUi d . ( 1 ) 

Here P k is polarization, u t is elastic displacement, f[ and are material 
tensors, (.. .),* = <9(.. .)/dxi , and Einstein summation rule is adopted. Note 
that usually strain tensor u t j = {u tJ +Uj t j)/2 is used instead of displacement 
gradients u h] (distortion tensor). But since is convolved with material 
tensors which are symmetrical in corresponding indices, it does not matter to 
write down u t ,j or Uij . Note also that here it is used definition of slightly 
different from that adopted in [TTj . pairs of indices ij and kl are swapped. 
This definition corresponds to [5) (9J ;TU| and more widely accepted. Certainly 
the opposite definition is also applicable. 

Thermodynamic potential of a body H is integral over body volume V 
from thermodynamic potential density. Using integration by parts its flexo- 
electric contribution can be represented in the form: 

a Six = J (■ Pk,i u i,j - p kUij,i) dV + j> d ijk UijP k dS , ( 2 ) 

where f k uj = is so-called flexocoupling tensor, dij k is surface 

piezoelectric tensor of special form: 



Hereinafter n; is unit vector normal to body surface S. Note that dij k is 
different in different points of the surface and it is symmetric in indices i, j. 

Thus flexoelectric part of thermodynamic potential density (jT]) can be 
reduced to the Lifshitz-type form 

T~Lfix ^ fklij (Pk,l'U , i,j P k Uij,l) (4) 

together with additional surface piezoelectricity described by piezoelectric 
tensor (J3]). Further effects of surface piezoelectricity are not described so 
that Lifshitz-type form (|4]) is used. 
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Certainly besides flexoelectric potential (j3J) one should take in account 
elastic energy. It will be clear from what follows that besides conventional 
elastic energy c ijk iUi d u k ,i /2 higher elasticity term v ijk i n mUi,j, n Uk,i,m /2 should 
be added. So that total volume density of the thermodynamic potential Hb 

is 

HB = 2 fa klUi * Uk ’ 1 + + fklij (Pk,lUi,j - PkUi^lj] • (5) 


Possible surface contribution to the thermodynamic potential is not consid¬ 
ered here. So that volume integral from ([5]) is total thermodynamic potential 
of a body. 

Starting from the equation ((5j) it is possible to write down the equations 
of elastic equilibrium and the boundary conditions for them. In fact, one just 
need to calculate the tensors 


1, p 

J-ij Qu- Cijkl'U'kJ, i 2 Jklij*k,l i 




( 6 ) 


@ijk 


dU 


B 


du 


i,j,k 


^ijnlk/nri .l.rn r , fnkij Pri • 


(7) 


and use the results of ca. In this way it turns out that “physical stress” ay,- 
which obey the equations of equilibrium 


cr, 




= 0 


( 8 ) 


are determined by equation: 


&ij Pij ®ijk,k Cijkl'U j k,l 'VijnlkmV j n,l,m,k T fklijPk,l • (9) 

Boundary conditions for (JHJ) can be written in different equivalent forms, one 
should only substitute (|6]) and (O to equations from [12]. Particularly they 
can be written in the form which coincides with that presented in [5]: 

®ijkBjTl k | g 0, (79) 

G’ij’B'j ®ijk,j^l k T @ijk,lBlTljTl k ®ijk^ljk\g 0 • (H) 

Naturally terms contained and Fi do not appear because surface contri¬ 
butions to thermodynamic potential is not described. 

It is obvious from (fTOj) and (J7|) that if Vijkinm —> 0 and fijkiPi\ s 7^ 0 then 
generally —>■ oo. So that the limit Ujjfcinm —7 0 is singular. This fact 
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explains the above statement that generally in the presence of flexoectric- 
ity one should take into account the higher elasticity (spatial dispersion of 
elastic moduli). Higher elasticity term can omitted only in very special (and 
physically doubtful) case when Pi — 0 on the surface of a body. In this 
special case one can use classical boundary conditions, additional terms in 
them disappear. 

Regarding the description of the direct flexoelectric effect, to obtain the 
corresponding differential equations and boundary conditions, one should 
add to (j4j) the following expression: 

Pp — ~^9ijklPi,jPk,l + 2°tjPtP] — EiPi , ( 12 ) 

where Ei is electric field. If surface contributions are not considered then it 
remains only to vary such thermodynamic potential in Pj and to obtain by 
standard way the following differential equations 


C^lkPl 9ijklPi,j,l fklijtkjjl Pfe 0 


and boundary conditions for them 


Qijkl Pi/J^l T 2 fklij'U'ijjW'l 


= 0 


(13) 


(14) 


It does not constitute a serious problem to take into account the surface 
contributions here. Usually such contributions contain only the polarization 
Pi but not polarization gradients P M -. This is why to find the change of the 
equation written above is a simple exercise in this case, in which we are not 
resting here. 


3 Exact solution for homogeneously polar¬ 
ized ball of isotropic dielectric 

When solving a particular problem it is convenient to use covariant formalism 
(see B2D- To transform equations from Cartesian form to covariant form one 
should only change partial derivatives to covariant derivatives and provide 
an index balance (only superscript can be convolved with subscript). So that 
equations of mechanical equilibrium and boundary conditions for them take 
the form: 

= ( 15 ) 
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cr“ 7 n. 


7 — 0 Q/37 ;) g?7. 7 


0 a/37 ngn 7 |^ ( = 0, 

+ <d a ^. s n 5 npn^ - 0 “ / 37 7 / 3 7 | = 0 


( 16 ) 

(17) 


Here n a is covariant representation of unit vector normal to body surface S, 
tensors a a/3 and 0"^ 7 are determined as follows: 


a af} = c a ^ s u r , s - v a ^u rM - £ + r 5af5 P r s , (18) 

0 a/37 = v a ^u eM ~ ^f 5iaP Ps ■ (19) 

It is also convenient to use curvilinear coordinate system of special class in 
which equation of body surface has the form x 3 = const. In this coordinate 
system tensor 7 °^ can be presented in terms of Cristoffel symbols |12j : 

7 / 3 7 = T £ l5 n 5 npn E + r^ s n s n e n y - T^ns . (20) 


Besides n a has the only component: 

n3 = 7F' 

For a case of a ball it is naturally to use spherical coordinate system: 

x — r sin 9 cos ip , 

< y — r sin 9 sin 'ip , (22) 

z = r cos 6 . 



Here x , y and 0 are Cartesian coordinates, ip = x 1 , 6 = x 2 , and r = x 3 are 
curvilinear ones. Equation of the ball surface has the form r = R, so that 
these curvilinear coordinates belong to the special class mentioned above. 

Instead of curvilinear indices 1,2, and 3 further it is used also the indices 
ip, 9, and r respectively. So the letters r , ip and 6 are excluded from notation 
of “running” indexes. Particular Cartesian indices are further denoted as x, 
y and z. 

Metric tensor components g a p , its determinant g , and Cristoffel symbols 
r “ 7 can be derived directly from ([ 22 ]) : 

2*2/) 2 1 4*2/i 

= r sin 9, g ee = r , g rr = 1 , g = r sin 9 , 


■ r*, = r*, = cot e, r * r = r* t =r-', r^ = -sm»cos». (23) 

. C = C = r" 1 , r^, = -r sin 2 0 , rj„ = -r. 
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Other components are zero. Since tensor g a p is diagonal, inverse tensor g°9 is 
diagonal also and corresponding components are simply equal to the inverse 
diagonal components of g a g . Unit vector normal to the surface has the only 
component n r — 1 . 

Further one should specify the form of the material tensors for isotropic 
media. It is well-known that elastic tensor of such a media is defined by two 
parameters, say its components c 12 and C 44 (Voigt notation is used), and can 
be expressed as follows: 


Cijkl Ci2&ij3kl V C 44 ( SikSjl V • (24) 

It is possible to show that flexocoupling tensor of isotropic media can be 
expressed analogously: 


fijkl f 12 ^ij^kl "f~ fAA^ik^jl "F $il$jk) • (25) 

For six-rank tensor of higher elastic moduli of isotropic media further it is 
used a simplified expression 

Vijlknm ^ 1 ( ^ij &lm $nk V 5ij5 n i5 m k V ^ij^nm^lk T /,- 6 j n Si rn V 

SilSj n Skm. V Si m 5j n Skl “I” ^in^jk^lm V Si n Sjl&km “I” Si n 8j m 8kl ) V (26) 

^2 (&ik$jl&nm. V 8{k8jm^nl “I” ^il^jk^nm *b 8^8jm^nk “I” 8j rn 8jk8nl V 
8im,8jl$nk) i 

where V\ and v 2 are constant parameters. 

Above the expressions of material tensors are presented in Cartesian com¬ 
ponents. To convert them to curvilinear components one should only change 
Kronecker delta to metric tensor with corresponding location (up or bottom) 
of indices. 

If the ball is polarized along c-axis with polarization P then polarization 
vector has two covariant components: 


P r = P cos 6 , 

(27) 

Pq = —rP sin 9 . 

(28) 


The equations above mathematically completely determine the problem 
under study. To solve this problem one must first separate the variables. 
Generally this is done by decomposition of the vector field u a in a series of 
spherical vectors. Note that since the covariant formalism is used here, spher¬ 
ical vectors do not coincide with their usual form [13]. However, the required 


form of these vectors can easily be obtained from conventional. Covariant 
spherical vectors can be defined as follows 



r 

1 lm\9 

\[K l 

+ i) 

v( 2 ) — 

i 


1 lm\9 

s/Ht 

+ i) 

y( 3 ) _ 
1 lm\r 

y 

1 lm i 


y(d . 

1 lm\r 

o, 

y ( 2 ) 

1 lui\r 




y V +1) 


imY ] m , 


-m 




Zm 1 


y ( 2 ) 

1 lull'll} 


/ 

, sin 6(Y lm ) t0 , 

V Z ( Z + !) (29) 


^(3) 

‘ lrn\Q 


y( 3 ) 


where Y) m is scalar spherical functions expressed in the usual manner in terms 
of associated Legendre polynomials. 

In the case considered here, i.e. when the polarization is homogeneous, 
expansion in spherical vectors is simplified radically. Indeed, it is easy to see 
that in such expansion of the polarization field P a there are only the terms 
with l — 1, m — 0 and top indices 1,3. Since the theory is linear, it follows 
that the expansion of the field u a contains only the similar terms. As a result, 
this leads to the fact that u r ~ cos 6 , ug ~ sin 9 , — 0. This is why one 

can use a substitution: 

u^ = 0, (30) 

u g = — r/ 2 (r)sind (31) 

u r — fi (r ) cos 9 . (32) 

Factor r added into (l3Tf to make radial function /^(u) analytical at r = 0. 
The consequence of using the covariant formalism is that u g is not a phys¬ 
ical displacement in the meridional direction, such physical displacement is 
u e\Jg e6 ■> while g ee = r~ 2 . Certainly the physical component of the displace¬ 
ment should not have a singularity at the ball center. Note that covariant 
form of Y^ e contains the same factor r. 

Thus, the problem is reduced to finding the radial functions f\(r) and 
/ 2 (r) . To find the corresponding differential equations and boundary con¬ 
ditions one should substitute in the above equations and make a 

very cumbersome but absolutely straightforward algebraic transformations. 
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As the result it yields the system of two differential equations 

v 3 (e 4 /r+4 eii - s efi + wo +- m)+ 

^(wrwwr-ewr+wo- 

v^f"' - 2f/" - 4^/2 + 12/2) = (33) 

c 44 £ 2 (2£ 2 /r + 4//( - 6/0 - C44 £ 2 (2£/' - 6 / 2 ) + 

Ci2^ 2 (e 2 /r + 2£/( - 2/0 - c 12 £ 2 (2£/' - 2/2), 

<w 4 /r+4/ 3 /2" - 4/ 2 /2 + 8/2)+u 3 (4/ 2 /r - 8/0- 

td (2 ef” - 4/2)+ v, (e 3 /r+ 4 / 2 // - 2^ - 4/0 = 

(34) 

c 44^ 2 (/^/2 + 2^/2 — 4 / 2 ) + C 44 / 2 (//( + 4/i) — 

Cl2^ 2 2/2 + Ci 2 / 2 (//( + 2/0 , 

and four boundary condition for it 

2^h (/0 + f '2 _ I8/2 + 34/2 + 8/f + 14/,' — 34/0 + 

+2u 2 (2/'" + 2/2 - 20/' + 36/ 2 + 4// + 16/( - 36/0- (35) 

-2 c 44 A 2 (/' - h + / i )| €=1 = fnR 2 P , 

R 2 [cu(2fi — 2/2) + (c 4 2 + 2 c 44 )/(] — 

U09/0 + 18/" - 46/; + 36.// - 14/'' + 32/' - 36/ 2 )- (36) 

2u 2 (3/(" + 6// - 22/( + 280 - 2/" + 16/' - 28/ 2 )| €=1 = P 2 /i 2 P, 

(18m + 12 v 2 )f? + 12m(3 ft - 4/! - 2/; + 4/ 2 )| ?=1 = 

(37) 

(/i2 + 2/ 44 )P 2 P, 
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2+ (/" + 2/2 — 6/2 + 2/,' + 6/1) + 


(38) 


4^ 2 (/" - 2/2 + 2/ 2 + 2/; - 2/01 ?=1 = / 44 P 2 P . 

Here £ = r/R is dimensionless radial coordinate, prime denotes the derivative 
in £, v 3 = (n 4 + 2 v 2 )R~ 2 , n 4 = (8u 4 + 4u 2 )P~ 2 . 

I 11 general, the system of two fourth-order equations (133|) — (T34|) requires 
eight boundary conditions, but here we have only four boundary conditions 
(l35|) — (1381) . However, as is usually in the case of such problems, the missing 
boundary conditions are replaced by the conditions that a solution should 
be analytic at zero. Therefore, the solution should be expressed in the form 
of series in non-negative powers of £. 

Practically it is more convenient to find the complete basis set of solutions 
expressed as such series, and to present a general solution as a linear combi¬ 
nation of such basis functions. If the number of basis functions is equal to the 
number of boundary conditions then the coefficients in the linear combination 
can be easily founded by solving a system of linear algebraic equations. 

In accordance with the above, we express a solution in the form: 

/.(e) = E«), ( 39 ) 

k =1 

00 

Bfci(0 = Z>Wnr- ( 4 °) 

n=0 

The constant coefficients a^n are found from the condition that Bki(£) obeys 
the system of differential equations (133]) - ([341) . This condition leads to the 
following linear algebraic equations for the coefficients (index k which is num¬ 
bering the solutions is omitted): 

[v;i(n 4 — 2n 3 — 9n 2 + lOn + 16) + 
u 4 (n 4 — 2 n 3 — 7 n 2 + 8 n + 12)]ai n 
+ [t> 3 ( 8 n 2 — 8 n — 16) — v 4 (2 n 3 — 8 n 2 + 2 n + 12 )]a 2n = 

[c 44 ( 2 n 2 - 6n - 2 ) + c 12 {n 2 - 3n)]ai n _ 2 - 
< [c 44 (2n - 10) + ci 2 (2n - 6 )]a 2n -2 , (41) 

[u 3 (4n 2 — 4n — 8 ) + n 4 (n 3 + n 2 — 4 n — 4)]a 4n + 

[u 3 (n 4 — 2 n 3 — 5 n 2 + 6 n + 8 ) — u 4 (2 n 2 — 2n — 4)]a 2n = 

. [c 44 (n + 2 ) + ci 2 n]ai n _ 2 + [c 44 (?r^ — 3 n — 2 ) — Ci 22 ]a 2n _ 2 • 
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From (HIT) one can see that the system of linear algebraic equations, even 
though it is infinite, has the characteristic structure: the coefficients with 
greater n are expressed in terms of the coefficients with smaller n. Therefore, 
one can fold ai n for any finite n step by step, one needs only to analyse the 
case of small n and define some of the cq n for these small n. 

Naturally for small n system of equations (FTTj) is degenerate, some aj n with 
small n can be defined arbitrarily. This arbitrariness leads to that it may be 
several functions Bki(£) and their normalization may be arbitrary. Detailed 
analysis shows that the system of equations (T4TT) is degenerate for n — 1,2,4 
only. It turns out five arbitrary constants, and so it is the same number of 
basis functions. However, one of this functions is physically meaningless and 
should be omitted (more on this see below). 

From an abstract point of view, arbitrary constants can be defined as 
anything, it is just needed that obtained basis functions be linearly indepen¬ 
dent. However, there is an important question: is it possible to choose these 
constants in such manner that the series representing at least some of the 
basic functions cut short, would constitute a finite sums? It turns out that 
this choice of constants is really possible. 

From the system of equations (14T| it is clear that in order to series cut 
short, it should be a 3 ; = a 4 * = 0. For this to be possible, in any case, 
right-hand sides of (1411) should be zero for n = 3,4. Together with (1411) for 
n = 0 , 1,2 these conditions lead to a system of linear algebraic equations 
which is quite amenable to analysis. The result is that there may be two 
basic functions which are finite sums. 

One of them is determined by the conditions a\o = a 2 o = const 7 ^ 0, 
a n >o i = 0. This function corresponds to the displacement of the ball as a 
whole. Such a displacement is physically meaningless and this basis func¬ 
tion should be omitted. Certainly such a displacement does not change the 
thermodynamic potential, so that this function indeed is a solution of the 
equations derived from a stationary condition of this potential. But it is 
meaningless solution. 

Another basis function, which is finite sum, is determined by condition 
(ci 2 —c 44 )a 2 2 = (3c 44 + 2ci 2 )ai 2 , other coefficients are zero. So that first basis 
function is determined as follows: 



(42) 
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other ai n are zero . 

It should be stressed that the function B\ obeys not only the equations 
(133|) — (1341) but also the equations of the classical theory of elasticity, which 
can be obtained from (1331) — (1341) by setting v 3 = v 4 = 0. It can be verified 
by direct substitution. 

Further analysis shows that there are three additional basic functions 
B 2 — B 4 which are expressed by infinite series and can be chosen, for exam¬ 
ple, as follows. 


' B 2 : 

_ 3C44 + 2Ci2 

012 — i , <222 —-012 , 

C44 — C12 

< (3-24 = 0 , 

(40u 3 + 60u 4 )ai4 = 

= (6C44 + 4 ci 2 )a -12 + (2C44 — 2 ci 2 )ci 22 , 
„ (3*0 (3*1 0 


coefficients with n 


coefficients with n 


> 5 are calculated by means of (1441) . 


B 3 : 

4u 3 + 3u 4 

<311 = 1 , <321 = ---—(3ll , 

4u 3 + 2v 4 

(3*0 = di 2 = (3*4 = 0 , 

= 3 and n > 5 are calculated by means of (TUI) . 


(43) 


(44) 


f B 4 : 

j _ 1 _ ^ V:i + /A Cl 

< d\ 4 — 1 , 024 — --(3'14 , (45) 

I V 4 - 4u 3 

( Ojo = Oil = Oi2 — Oj 3 = 0 , 

coefficients with n > 5 are calculated by means of OH]). 

Graphs of the basis functions Bik defined above for some set of parameters 
are shown in FigjTl Turn our attention that in contrast with Bn the functions 
B 2% — B 4l are concentrated near the surface and decay rapidly in the ball 
volume. 

Above there is explicit representation of four basic functions B\ — B 4 . 
It only remains to substitute (j39l) in the boundary conditions (j35l) — (j38]h 
and from the resulting system of four linear algebraic equations to find nu¬ 
merically Ck ■ Thus, the problem is solved. Result of numerical calculations 
for some parameters is presented in the Fig]2] 
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Figure 1: Basis functions for B — 1-10 5 , P = 1, c 44 = 1.1 • 10 12 , c 12 = 
3.4 • 10 12 , / 44 = /12 = 1 • 10- 3 , V! = 0.2, w 2 = 0.1. 


Note that the graphs shows that, except for a thin layer near £ = 1, 
the curves are parabolic, the more so the less V\ and u 2 are. In the bulk 
of the ball the only function £hi(0 remains with a good approximation. 
Thus, for a small Uj, except for a thin surface layer, solution approximately 
obeys the equations of classical elasticity theory (see discussion of function B\ 
properties above). This observation is used in the next section to construct an 
approximate method of flexoelectric deformations calculations. This method 
is needed to describe more complicated geometry then spherical geometry of 
a ball. 
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Figure 2: Result of the calculation for R = 1 • ICR 5 , P — 1 , C44 — 1.1 • 10 12 , 
C 12 = 3.4-10 12 , f 44 = f \2 = T10 -3 . Solid line is for V\ = 2-ICE 1 , V 2 = THE 1 , 
dashed line is for v± = 2 ■ ICE 2 , V 2 — 1 • ICE 2 . 

4 Approximate method to calculate flexo- 
electric deformations of finite size bodies 

Exact (within the framework of a continuum media theory) equations of 
elastic equilibrium used in previous section are partial differential equations 
of the fourth order. They are too complex to be used in the cases interesting 
for applications. Even for a homogeneously polarized ball of isotropic media 
solution of them is cumbersome and requires a functions expressed by a series 
in powers of radial coordinate. So that to find flexoelectric deformation in 
more complicated cases one needs an approximate method. Such a method 
was proposed in HD] and it is described in this section. 

When solving the problem for a ball (see previous section) we have seen 
that the solution is a sum of two parts. The first part corresponded basis 
function B\ obeys the equations of the classical theory of elasticity, the second 
(contribution of remaining basis functions) is a non-classical part which, for 
small v , is concentrated near the surface and rapidly decays inside the 
body. Hereafter the first part is called a volume or a classical part, and the 
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second one — a non-classical or a surface part. 

It is naturally to assume that this property should be preserved in more 
general cases. The aim is to use this observation for construction of approx¬ 
imation that is useful for solution of more general problems than the case of 

a ball. 

By means of (UHl) equations (ITS]) can be represent in following form 
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u rAl , + P^P-rM - = 0 . 


(46) 


Boundary conditions for this differential equations here are used in the form 
different from used in previous section. Particularly they are 


0" 33 = 0 , (47) 

cr"3 - 0^ 3 ;/J + 0“ 33 3 + ©“(^r 3 ^ = 0 . (48) 

ffere the indices enclosed in parentheses runs only the values 1 and 2 , if there 
are multiple indexes in one pair of parentheses then these indices are not equal 
to 3 simultaneously. Equivalence of this form of boundary conditions to (fT71) 
is proved in [12]. Certainly it is assumed coordinate system where equation 
of the body surface has a form x 3 = x|, is some constant. 

According to the observation mentioned above the solution of the equa¬ 
tion (1461) is presented in the form w 7 = w 7 + w 7 , where n 7 obey the classical 
equations: 

C a ^ 5 U r 8;P + r Sa PP r ,8;(3 = 0. (49) 

From (ITS]) and (ITS]) it follows that h 7 obey differential equations: 


c“^ s u rS .f - v al> ~' l ‘ ( u rMil 0 = r j. ( . S 0 . 


(50) 


Substitution w 7 = n 7 + h 7 and some transformations also yield that the 
boundary conditions (TT71) and (ITS!) take the form: 


v 


a3eS3£ 


_ * c63a3 p 


a3e53(~ 
v U e-M, 


(51) 
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e” 37i fi 7 ;J + r la3 P r ,5 - ~ U, MiW + 


v°W‘^u eM . m + v“V>‘Wu cM r 3 m7) = 0 . 


(52) 


Equations above are exact. In order to turn into approximation one 
should note that real values of are small. It is convenient to assume 

that all of them are proportional to a scalar v —> 0. It is assumed also that 
all c a ^ 5 is proportional to a scalar c. If v —> 0 then the right-hand side of 
dSH]) and the second term on the right hand side of (T5TT) can be neglected. 
Indeed, since ft 7 obey the classical equation (149|) . in the case of sufficiently 
smooth surface and in the absence of polarization gradients tend to infinity 
in the limit v —> 0 , the classical part can not have large gradients which 
can compensate smallness of v " /37<5e< ’. As for the non-classical part w 7 , the 
situation is different: here such compensation is possible, but only if the 
derivatives are in x 3 . Moreover, in a thin layer near the surface covariant 
derivatives in x 3 can be replaced by the usual derivatives, T-terms give only 
small corrections here. Thus, (I50|) and (l5Tf can be approximately replaced 
by 


C“ W « 7 ,3,3 - «“ 37333 * 7 ,3,3,3.3 = 0, 


(53) 


«3e333 - 1 r63a3 

V U £ , 3,3 — 2^ 


Pa. 


(54) 


Naturally decomposition u 7 = -u 7 + 'ii 7 is not unique. By requiring that h 7 
obey (H§p . this non-uniqueness is restricted, but is not eliminated completely. 
It is clear that h 7 does not contain a non-classical part, but it does not 
mean that u 7 does not contain a classical part. To eliminate remaining non¬ 
uniqueness one should require that u 7 exponentially decay inside the body. 
Under such condition u 7 obey not only (1531) but also 


„o373 


U 7 — V 


q37333 


^ 7 , 3,3 — 0 . 


(55) 


Equations (l55|) are a system of three ordinal differential equations, the 
dependence on x l and x 2 is parametric here. Moreover, in a thin layer near 
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the body surface one can assume that the coefficients do not depend on x 3 , 
they are approximately equal to the surface values. Solution of such a system 
is easy. In a standard way one should find a fundamental basis set of solutions 
in the form w 7 = u y e x ( x ~ x s '>. Obviously equations for amplitudes fi 7 are 

c“ 373 u 7 = A V 37333 u 7 . (56) 

Equation (1551) is a standard generalized eigenvalue problem for symmet¬ 
ric positive de fini te 3x3 matrices. Therefore X n and u" are calculable by 
standard way. If the body bulk corresponds to x 3 < x 3 s for definiteness, 
then A n equals the positive square root of the n-th eigenvalue, and general 
representation of non-classical part is 


fi 7 = E a n u: f e x ^- x s ), 

n =1 


(57) 


where only the coefficients a n are unknown. The latter can be found easily 
by (154)) which yields a simple system of linear algebraic equations: 


E A^ 33 Xa„ = i f‘**p s . (58) 

n= 1 ^ 


Note that the smallness of u" 37333 in (1581) is compensated by A^, and 
using (156|) this system of equations can be rewritten as follows: 


£ c^u^an = \f 53a3 P s . (59) 

71=1 ^ 


Thus, the non-classical part of clastic displacements is found completely in 
explicit form. 

As it follows from above, to find completely the non-classical part of 
clastic displacements one needs only the boundary conditions (1511) . The 
remaining boundary conditions (1521) yield the boundary conditions for the 
classical equations ()49l) in this way, one should only substitute known u 7 to 
(1521) and select the terms of the appropriate order of accuracy. 

However, within this derivation of the boundary conditions for the equa¬ 
tions (149|) a singularity arises requiring additional analysis. This singularity 
is related to the fact that it appears a subtraction of two close but large 
terms in which it is impossible to use the approximate expression for u 1 de¬ 
rived above. In these (and only these) terms it is necessary to consider a 
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corrections to such h 7 . The appropriate analysis is made below and while 
w 7 means the exact non-classical part of elastic displacements rather then 
approximate one defined by ([57]) . 

According to the above, all the terms in (|52j) . where gradients of u 1 are 
convolved with , should be omitted for the same reasons as in the case 

of the equation (f5U[) . Further it is convenient to introduce the notation: 


s“ = -c a3 ^u r ,s + v a3 ^u rM , £ + v a ^ 5 ^u rMm - 


v a(P)e6( t)C - p3 
V U e-M L 03 )( 7 ) 


(60) 


With this notation the boundary conditions to the classical equations can be 
written as follows: 


C a ^ 5 u r ,s = s a - r Sa3 P r ,s -\f 53a{m Psm + lf 5 h)a{ 0 ) P s r 3 mi) . (61) 


Now all the terms to be further simplified are contained in s a . Carrying 
out simplification, first of all one should leave in the last term only the partial 
(not covariant) derivatives in x 3 . The reason for this is the same as above. 
Next, one should express the remaining covariant derivatives in terms of 
partial derivatives and T-terms. The exact expression for the third covariant 
derivative of the vector is very cumbersome. But keeping in mind that these 
third covariant derivatives are convolved with v ^ one can om q all the 
terms where initial vector is not differentiated at least two times. So it turns 
out: 

^ a ;/5;7iC ~ 7,C — U 8,P£^'a'y — U 6,'Y,C^'aP ~~ a(~ 

(62) 

h a ,5,(P — . 
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With this equation and simplifications mentioned above we get: 

q" — —r a3 ' yS 'Ti r 4- . _ ..o&ySeCf pP _ 

b C -(5 nr U ^7,(5,^,6: " 7^ 

oiSrjSeiif pP — o, a3 'Y S <f, c ,pP _ pP _ 

“ “p,C,A 76 “ U P,S, C 1 7£ “ U l,P,e L 5C 


“ ' ’“p,^ 1 75 - “ ' “w 1 * - 

F 3lfeC «„, c r; e - F^fl 1jV q e + </-«w<4 lAC0S) - 


v a (Ph S3 Cf, pP _ v a (Ph 63 C 7 7 A pP _ v a(Phfi 3 Cf, pP . 
V U 1 ,P,(P) L <5C 1 “TSPiC 1 <5(/3) 1 WAp 1 C(/3) 


v a(^) £ 3(7)3 4)3j3r 3^ w 

Again, using the fact that to compensate smallness of ^“Ar&C one needs at 
least two derivatives in x 3 , this equation is further simplified: 

s a = _ c ^u rS + v -^u lM e - V^^Up^T^- 

n a 37 < 533 Mp ) 3 i 3 r^ - 7 , a 373 e 3 u P) 3 ) 3 r$ e - n Q 3753 ^ 7 , 3 , 3 r 3 c - 

a3 7 &3 7 p3 a3 7 3£C 7 F 3 _j_ „,a 03)7333 7 ^ \ 

^ “ 7 , 3 , 3 ! fe~“ s “ 7j3)3 i f£ + V “ 7 , 3 , 3 ,(/?) (64) 

7 ,a(£b 333 7 7 pP _ a(/3)7<533 ~ „p 3 _ 

v “p,3,3l 7 (yj) “ “7,3, 3i <5(/3) 

a(/3) 7 33C 7 p3 a(/ 8 )e 3 ( 7)3 7 p3 

“ “ 7 , 3 , 3 ! c(/?) - V “e,3,3! ( 0 X 7 ) • 

It is convenient to introduce yet another notation: 

^op _ y o3733CpP^ _|_ v a37<533pp^ _|_ v a3 7 3£3pp _|_ v a3p53Cp3^_|_ 

^a 3 p 5 £ 3 p 3 £ + v a 3 p 3 £Cp 3 £ + v a(/ 3 ) 7 333 pP^ + ( 65 ) 

v a(/3)p533p3^ ) + ^(flpSSCpS^ + v a(/3)p3 ( 7)3p3^^ _ 

Using the symmetry properties of it can be rewritten as follows: 

/i a/3 = (u a373,53 + 2n a37633 )r^ + ( 2 v a3 ^ 53 + n Q3/373,5 )r 3 5 + 

(66) 

2 v a3/37(5)3p3 (5) + ^373(5)3^ + ^(^S^SpS^ _ 
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In this notation the expression for s a can be written as follows: 


= -C a3l5 U lt s + C a3 ^U £ T £ l5 + V a3 ^U^ M ,e+ 
v a ^ 333 u Mi 0 ) -h a %, 3 , 3 . 


(67) 


Now one can proceed to the analysis of the singularity mentioned above. 
It appears in the first and third term of (l67j) . when all the derivatives are in 
x 3 . For all other terms there is no singularity and one can use the expression 
(1571) for u 7 . It is also clear that singular terms exactly cancel each other out 
because if approximate expression (1571) is used. This is due to the explicit 
form of this expression. One might think that these terms should simply be 
excluded, but in reality the situation is more complicated. 

The fact is that —c" 373 'U 7)3i3 + v a3l333 u lt3 ^ 3t3 is equal to zero only for the 
main approximation for u 7 . But the corrections to this approximation may 
result in the fact that this expression becomes finite. Henceforth we denote 
such corrections as w-y keeping notation u 7 only for the main approximation. 
In this notation it turns out: 


s“ = -c “ 373 ™ 7 , 3 + v a3 ^ 333 w M3 - c a3 ^u^ {S) + c a3 ^u £ T £ lS + 
V a3 ^% (M ,e) + V a ^ 333 U^ 3, 3 ,( 0 ) - h a %, 3,3 • 


( 68 ) 


To compensate the smallness of v a/3lS£ ^ in the fifth term of this equation there 
should be at least two differentiations in x 3 . Thus, after some reductions one 
can rewrite (168|) as follows: 


s“ = -c “ 373 «; 7 , 3 + w a 37333 w 7i3 , 3 , 3 - c a3 ^u lXS ) + c a3 ^ 5 u £ r £ lS + 

2 ( l , a3 ^ 33 + U q 373 ^ 3 K 3 , 3 , (/3 ) - h a %, 3 , 3 . 


(69) 


It is clear that in (RTJj) one should use only corrections w 1 of order v 1 ^ 2 . 
Only such corrections yield a contribution of the same order as the other 
terms. To find the corresponding contributions, we should keep in the exact 
equation 


c a ^ 6 (u. 


+ w 7 ) 


5;/3 


V 


aP'ySeC, 


Ury + ia 7 ) ; 5;C; e;i a 






(70) 


only the terms of the order v 1 / 2 . Indeed if w 1 ~ u 1//2 then the main terms 
w 7j3i 3 and v a37333 w 7 , 3 , 3 , 3,3 have just such an order. Note that although there 
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are terms rt 7)3)3 and v a3 ' y333 u lt3 ^, 3,3 in the equation, which are of the order 
~ v ^ 1 , now they cancel each other completely due to the fact that now w 7 
denotes the main approximation. Right hand side of (1701) is ~ v , so it should 
be omitted. The terms c^^w^s-p) and v a ^ Se ^w 7 (.s-Q S -p) should be omitted 
by the same reasons. Thus (1701) is reduced to 


- C“ 31 V,3.3 + <'“ 373S X 3 ,3,3,3 = C“' 53i fi 7 {m ~ . 


(71) 


Equation (ITT)) requires further simplification: in the right-hand side we 
should keep only the terms ~ n -1 / 2 . For this in the Erst term there should 
be one partial derivative in x 3 , and in the second term there should be 
three derivatives in x 3 . This is why we simplify second and fourth covariant 
derivatives as follows: 


h 7; 5 ; /3 ~ "h 7 ,<5 1( 8 Ur y,p^SP ^P><5^ 7 /3 ^'Pt^'yS i 


(72) 


^7;i5;C;£;/3 ~ ^7,S,C,e,p ~~ ^p,S,C,e^jp ~~ ^"y,p,C,^sp~ 
V"V,8,P,e^'cP ~ ^yACjP^e/S — ^p,5,e ,^ 7 £ — ^ p,C,e,P^8~ 
UpM^e ~ ~ ^7,/5>C,£^& ~ ^"1 Ap,P^(,£ ' 

Eventually the simplified equation (1711) takes the form: 

-C“ 3 " 3 U; 7 ,3,3 + ^“ 37333 ^ 7) 3,3,3,3 = (c^ + C^)u^- 
c a ^ s Y]pU^ - (c“^ 3 + c“ 37/3 )T^h £ , 3 - 2 (u q37(/3)33 + 
v a ^ 333 )u , m ^ 3 + 2(u^ 7333 r + u“ 37533 r^)h Pi3 ,3,3+ 

(4u ^ 7 i33 r 3^ + ^3*3^ + ^53^3^^ 


(73) 


(74) 


Taking into account that near the body surface the Christoffel symbols 
and material tensors can be considered as constants, the equation (1741) has 
a very specific form: all of its terms are differentiated with respect to x 3 
at least once. Therefore, there exists a particular solution of this equation 
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which obeys an equation with one less derivation: 


-C a 3 ^ 3 W y , 3 + V a 3 ^ 333 W lt 3 , 3,3 = (c° 37(/3) + C a ^ 3 )u im - 


c Q/? 7 < 5 r^u 7 - (c a/?73 + c q 37 / 3 )T^« £ - 2 (u q37 ^ 33 + 

(75) 

v a ^ 333 )u im 3,3 + 2(u^ 7333 r^ + n Q37533 r^)n p ,3,3+ 

(4 v a ^ 7 ^ 33 r 3 ^ + v «^ 3 r^ + n a3 ^r 3 ^)ti 7 , 3) 3. 

Note that in the left-hand side of (175]) there is exactly the expression that 
is needed in (169]) . The general solution of the inhomogeneous differential 
equation is the sum of a particular solution of the inhomogeneous equation 
and the general solution of the homogeneous equation. But the remarkable 
fact is that any solution of the homogeneous equation decaying away from 
the surface when substituted into —c a3 ' y3 w ^ /t 3 + r’“ 37333 iu 7 , 3 , 3, 3 yields zero. So 
it is quite enough to consider only a particular solution, and one can simply 
replace — c a3l3 w 1)3 + n Q 37333 ta 7 , 3 , 3, 3 in (1551) by right-hand side of (1751) . ffaving 
made the replacement, we obtain s a in the form: 


s“ = -c a37(5) u 


7,(5) + C a3l& U 


e^S + 


C 


_|_ c a(/?b3 


) U 7 XP)~ 


c a ^ 5 r ^ 7 - (c a/?73 + c a3 ^)Y £ lfi u £ - 2(v a3 ^ )33 + 

„, a ( P ) 7333\ - 1 a37(5)33 7 , 9 a373(/3)3 7 

V )«7,(/9),3,3 + 2n M 7 ,3,3,(5) + 2u M 7 ,3,3,(/3)- 


(76) 


^“0,3,3 • 

Elementary transformations lead (1761) to a rather simple form: 

= c a ^ 3 u im - c a ^ s r 3 ^ 7 - c“^ 73 r^« e - /l °%, 3i3 , (77) 

where /r a/3 is redefined as 


/z «/3 = u a(7)^3(i)3 r 3^^ _ ^a3 7 333 r ^ _ _ ^ 363 ^ _ ( 7 g) 


This equations together with (16T1) completely solves the problem considered 
in this section. 

Remember that in used coordinate system unit vector normal to the body 
surface n a has the only component n 3 = 1 /. So that with respect to the 
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classical part of strain, c a3 ' /S u Tt s actually defines the external force density 
a al3 np on the body surface, one should only multiply it by 1 /\fg iA - This is 
why differential equations (149|) appended by (l6Tf correspond to a standard 
problem of classical theory of elasticity for a body under external forces on 
its surface. Methods of solution of such problems are well-known and do not 
require special consideration. 

It should be emphasized that physically there are no external forces on 
the surface. Forces mentioned above are purely formal in nature and describe 
the interaction of non-classical and classical parts of deformation. 


5 Comparison of exact and approximate so¬ 
lutions for a homogeneously polarized ball 


It worth to find flexoelectric deformations of homogeneous polarized ball by 
approximate method described in previous section and compare this defor¬ 
mations with ones known from exact solution. This is a test of approximate 
method at least. 

For isotropic material the matrices c“ 373 and v a37333 are diagonal. So that 
in order to find u 7 one does not even need to solve the generalized eigenvalue 
problem (l56]h the result is obtained at once: 


where 


- ./12 + 2/44 X r R (£- 1) 

u r = —7 - -PcosOe r ^ ’, 

2( c 12 + 2 c 44 ) 


Ug = 


-R—P sin 9e Xem - 1] , 


2 c, 


44 


_ fly 0 7 


A g — 


\ r — 


C44 


Vi + 2 v 2 ’ 


I C 12 + 2c, 


44 


(79) 


(80) 


(81) 


9rq + 6 v 2 

For comparison with the exact solution it is useful to present in terms of 
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the functions f x and / 2 : 


u r 

cos 6 


/12 + 2/44 pc \ r R ( £ - l ) 
2 (ci 2 + 2C44) 


(82) 


/ 2 = 


u e 


1 

1 


fu 


p^eR^-l) 


(83) 


r sin 6 £ 2 C 44 

Next step is to find a a3 n 3 = c ol3lS u~ rt sn 3 on the ball surface. To do this, 
one should just use the equations of the previous section. It turns out 


a 33 n. 


P cos 6 

2 C 4 4 (/i 2 + 2 / 44 ) nf 

S R 

1 O V 44 

c 12 + 2 C 44 


(84) 


P sin 9 

/12 

c 12(/l2 + 2/44) 

s i ? 2 

2 

2(ci2 + 2C44) 


(85) 


Certainly d 13 = 0. This is obviously from symmetry but can be also obtained 
by direct calculations. 

To solve the differential equations f)49l) actually is not necessary here. It 
is clear that in the terms of functions f\ = u r / cos 9 and / 2 = — uq/(t sin 9) 
the solution is proportional to the previously defined function B x . Eventually 
the validity of this statement can be verified by direct substitution. So that 
one just needs to find the coefficient of proportionality, it is easily done by 
using (T84l) or (J85J). One can use any of these conditions, they both give the 
same result: 

7 _ PjCnlAA ~ O4./12) (02 ~ C44) ^ 2 /ggx 

£44(02 + 2C44) ( 3 cq 2 + 2C44) 


7 _ P( c 12/44 — £44/12) ( 2 ci 2 + 3c 44 ) 2 /gyx 

044(012 + 2C44) ( 3 Ci 2 + 2C44) 

It is also useful to derive the expression for the classical part of the elastic 
displacement in Cartesian components as a function of the Cartesian coordi¬ 
nates. Direct conversions yield 

( u z = aiz 2 + a 2 (x 2 + y 2 ), 

< u y = (a x - a 2 )yz , (88) 

[ u x = (ai - a 2 )xz , 
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where 


(89) 


'.'i = 


P\Cl 2 fu ~ C44/12) (Cj2 — C44) 

-R 2 C 44 (Ci2 + 2C44) ( 3 Ci 2 + 2C44) 


02 — 


-P(ci2/44 — C 44 / 12 ) (2Ci2 + 3C44) 
i? 2 044(012 + 2C44) (3Ci2 + 2C44) 


(90) 


Now we can find /, = /) + f) and compare it with the results of ex¬ 
act calculations. This comparison is shown in Figj3j This figure shows a 
good agreement between the approximate and the exact solution, the smaller 
higher elastic moduli match those better. It should also be emphasized that 
the latest version, when the difference is badly distinguishable, corresponds 
to most physically reasonable values of higher elastic moduli . Thus, the 
approximate method works very well at least for a homogeneously polarized 
ball with the given parameters. 
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Figure 3: Comparison of exact and approximate solutions for a ball. R = 
1 • 10~ 5 , P = 1, c 44 = 1.1 • 10 12 , C 12 = 3.4 • 10 12 , f u = fu = 1 • 10” 3 . a - for 
V\ = 2.0 • 10" 1 , v 2 = 1.0 • 10" 1 ; b - for Vl = 2.0 • 10" 2 , v 2 = 1.0 • 10" 2 ; c - for 
Vi = 2.0 • 10” 3 , v 2 = 1.0 • 10~ 3 . Solid line is approximate solution, dashed 
line is exact solution. 






6 Flexoelectric bending of homogeneously 
polarized circular rod 


In this section we apply the above-described method to find the bending of 
a homogeneously polarized circular rod of an isotropic material [2j. Such 
a rod may be polarized in different directions. However, due to the fact 
that the flexoelectric effect is linear, we can restrict ourselves to the case of 
longitudinal and transverse polarization only. The response of the rod to the 
polarization of an arbitrary direction obviously is a superposition of responses 
to the longitudinal and transverse polarization. Besides, axial symmetry 
allows us not distinguish different directions of the transverse polarization. 

As stated above, the classical part of strain is determined by the formal 
forces on the body surface. The surface of the rod consists of a cylindrical 
surface and end surfaces. First we consider the cylindrical surface of the 
rod. Here it is natural to use a cylindrical coordinate system. Let cylinder 
be located along the longitudinal axis corresponded to coordinate x. This 
coordinate is also denoted as x l . The other two curvilinear coordinates are 
the angular x 2 = 0 and radial x 3 = r. Angle 0 is measured in respect to the 
direction of the Cartesian z-axis. Equation of the surface is x 3 = R where 
R is rod radius. So that this coordinate system belongs to the class that is 
needed in accordance with the method of calculations. 

For such coordinates simple geometrical reasoning lead to the following 
equations defining the relationship between Cartesian and curvilinear coor¬ 
dinates: 

i 

x — x , 

< y — r sin0 = x 3 sin a; 2 , (91) 


z = r cos 0 = x 3 cos x 2 . 

By means of simple differentiation it is easy to find the metric tensor, its 
determinants and the Christoffel symbols: 

( 9ll 9xx 1 j 922 9w T i 933 9rr 1 j 9 T i 



~P<A _ ~P<A _ 
1 r(j> L 4>r 


l 

r 


(92) 


The remaining components are zero. Since the metric tensor is diagonal 
to find the inverse metric tensor g al3 is easy: one just need to take the 




reciprocals of the diagonal components. Note also that the vector normal 
to surface has the only component n r = 1 . 

By making perfectly straightforward calculations similar to the case of a 
ball described above, the following expression can be obtained for the non- 
classical part of the elastic displacement u a near the cylindrical surface. If 
the polarization P is longitudinal along the x-axis then the expression is as 
follows: 


U x U\ -P 


/44 




2 c. 


44 


U ( j > =U2= 0 , 


(93) 


U r =U 3 = 0 . 


For transverse polarization along z-axis u a is determined by the following 
equations: 

u x 0 , 


U<f, = 


-PR— sm0e A * fl(€-1) , 
2 c 44 v 


D f 12 + 2 / 44 X r R(£—l) 

u r — P —7 - -cos0e r ’ 

2 (c 12 + 2 c 44 ) 

Here ^ — r/R is dimensionless radial coordinate, 


'Vt, Af/, 


C 44 


ui + 2 u 2 ’ 


(94) 


(95) 


At' - 


I C\2 ~\~ 2C44 


(96) 


9u 4 + 6n 2 

Further, direct application of the equations of section H] gives the following 
boundary conditions for the classical part of the displacement. For transverse 
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polarization they are 


d 13 n 3 


0 , 


d 23 n 3 


d 33 n 3 


P sin ( 

R 2 


P COS( 

R 


C12CA2 + 2/44) 
C12 + 2C44 


fl2 


(97) 


c 44(/l2 + 2/44) 


C 12 + 2 a 


— fu 


44 


If polarization is longitudinal then all (T q3 are zero. 

Boundary conditions (197|) can be also expressed in terms of the Cartesian 
tensor component. Elementary transformations yield: 


Oils'll5 — 0 , 




■ _ -P(c44/i2 — C12/44) 

“ i?(c 12 + 2 c 44 ) 


• sin 20, 


(98) 


_ -P(C4 4 /l2 ~ £12/44) 

-R(Cl 2 + 2 C 44 ) 


cos 20. 


In this form of boundary conditions it is immediately clear that the total 
force acting on the cylindrical surface is equal to zero, the integration over 
the angle turns these expressions to zero. It also can be easily calculated that 
the total bending moment of these forces vanishes while integrating over the 
angle 0. 

Thus, not only the sum of formal surface forces, but also the sum of the 
bending moments of these forces, acting on a small part of the cylindrical 
surface length, are equal to zero. This means that the bending moment does 
not change along the homogeneously polarized rod. We emphasize that even 
so, the bending moment appears, but entirely by the boundary effects at 
the ends of the rod. So that in the case of homogeneously polarized rod 
the problem is reduced to the standard problem of the classical theory of 
elasticity [13], i.e. to the determination of rod bending under the bending 
moments applied to its ends. 

Certainly the surface forces (197|) slightly deform cross-section of the rod. 
Since the main effect is the bending of the rod, we do not discuss this defor¬ 
mation. Note only that qualitatively it is similar to the deformation of the 
ball in the meridian cross-section. 
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To calculate the bending moment appeared at the ends of the rod, we 
need to specify the shape of these ends. We emphasize that we can not limit 
the rod by planes because this gives sharp edges while the theory requires 
smooth surface. We should smooth out these sharp edges, say by a quarter 
of toroidal surface, or assume that the rod terminates by halves of the ball. 
Calculations were made for both these variants, and it was obtained the 
same result. Omitting the details of these rather simple calculations, we 
note only that in the case of spherical ends one can use the results of the 
previous section. For edges slightly smoothed by the toroidal surface one can 
approximately describe a small part of this toroidal surface as cylindrical one 
and use the equations given above in this section. Eventually it turns out 
that the bending moment is 

M _ P7ri? 2 (c 44 /i 2 - 012/44) 
v ~ C12 + 2 c 44 

The other components are equal to zero. It is implied here that the rod is 
transversely polarized along £-axis and the equation is written for the end 
in the positive direction of x-axis. For the other end of the rod M y has the 
opposite sign. 

As is already clear from symmetry, for the longitudinal polarized rod there 
is no bending moment at all, this can be also proved by direct calculations. 
Such calculations also show that the total force acting on each end of the rod 
is zero as for longitudinal and for transverse polarization. 

Equation (l99|) actually solves the problem of flexoelectric bending of the 
homogeneously polarized rod. One should only substitute this bending mo¬ 
ment into standard equations from textbooks (see [12] for instance). 

7 Flexoelectric bending of a homogeneously 
polarized circular plate 

In this section we apply the above-described method of approximate calcu¬ 
lation to find the flexoelectric bending of a thin circular plate with a radius 
R and thickness h, uniformly polarized normal to its plane p3]j. Material of 
the plate is assumed isotropic. For definiteness, we assume that the average 
surface of the plate lies in the coordinate plane OXY. 

Calculations for the plate are generally similar to those made in the pre¬ 
vious section for the rod. However, most part of the plate surface is flat, so 
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that one can immediately conclude, without any calculations, that the formal 
forces on this part of the surface are zero. As in the case of the rod, bending 
of the plate is determined by edge effects which are discussed below. 

Thus, for most part of the plate, except for the edges, the classical part of 
the plate deformation is determined by standard equations without surface 
loading. It is well known [Uj that under such conditions the displacements 
of the average surface of the plate ((x,y) obey two-dimensional differential 
equation 

AAC = 0, (100) 

where A is two-dimensional Laplace operator. The components of the strain 
tensor can be expressed in terms of ( as follows: 


^xx %C,x,x ; ^yy %C,y,y i ^xy %C,x,y i 

Hxz — ^yz 0 , Uzz % (C,x,x ~t~ C,y,y') 


( 101 ) 


C \2 T 2c, 


•44 


Using polar coordinates it is easy to find an axially symmetric, regular at 
origin, solution of the equation (HOOD : 


. Gr 2 G. 9 9 . . 

C {x,y) = -— = ~{x + y ). (102) 

This solution contains a single integration constant G which is nothing but 
the curvature of the plate. 

To find integration constant G one needs first calculate the components 
of the strain tensor by means of (11011) and second calculate the components 
of the strain tensor by means of standard equations a t j = CijkiUij . Thereafter 
it becomes obvious that there is the linear density of the bending moment M 
on the boundary of the plate. For instance, at the point of intersection of the 
plate boundary with the coordinate axis OX this bending moment density 
is 

+h/2 

1 ,/f / ( \ ,1 h C 44 (3ci2 + 2 c 44 ) 

M = / a xx (z)zdz = — G -—-. (103) 

J D C 4 2 + 2 c 44 

—h/2 

Due to axial symmetry M is the same at all other points of boundary. 

On the other hand M can be calculated in terms of formal forces arising 
due to flexoelectricity on curved surfaces on the plate boundary. As in the 
case of the rod, it is necessary to smooth out sharp edges using one-quarter 
of toroidal surface (on each sharp edge) or one-half of the toroidal surface for 
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the whole boundary. A small part of a toroidal surface can be considered as 
a cylindrical one, so that we can apply the equations (TfJST) (coordinate system 
should be rotated). For both variants of sharp edges smoothing it turns out 

M = F/ ^/l2-Cl2/44) ( 104 ) 

C 12 + 2C44 

It only remains to equate the two expressions for M and express the curvature 
of the plate in terms of polarization P : 


q _ 6-P(c44/i2 — C 12 / 44 ) 
h 2 C44(3Ci2 + 2 C 44 ) 


(105) 


It is interesting to compare the equation (I105p with the result obtained in 
[4] by direct minimization of the plate energy. Elementary transformations 
of the equations from [4j yield the equation different from (1 105 11 only in that 
there is 12 instead 6 in the numerator. However, it should be kept in mind 
that in [4] it is considered the case when in a thin layer near the plate surface 
polarization falls to zero. If, in accordance with calculations presented here, 
we modify the calculations [0] to the case, when the polarization is strictly 
homogeneous, then it appears exactly the equation (11051b 

A similar dependence on the boundary conditions imposed on the po¬ 
larization also holds for direct flexoelectric effect, it is discussed in detail in 
another paper [T6] , 


8 Conclusion 

Above we have considered the continuum theory of the flexoelectric effect in 
the finite-size bodies. In the part of description of the converse flexoelectric 
effect, this theory turns out to be quite complicated. This is due to the fact 
that the independent variables are the elastic displacement. Flexoelectric 
part of the thermodynamic potential depends on the second spatial deriva¬ 
tives of the elastic displacement. Moreover in the case of flexoelectricity the 
elastic energy should be considered taking into account its dependence on the 
second spatial derivatives of elastic displacement (elastic spatial dispersion), 
or the theory is self-contradictory in the general case. All of this leads to 
several theoretical problems. 

The Erst theoretical problem is that derivation of the boundary condi¬ 
tions for the differential equations of elastic equilibrium becomes non-trivial 
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in such a case. While the volume integration by parts, which is necessary to 
express the variation of the thermodynamic potential in terms of independent 
variates, its appear the surface integrals containing the gradients of these in¬ 
dependent variates. Needed transformation of these surface integrals require 
special mathematical tools. For this purpose it was offered to use additional 
surface curvilinear coordinate system or to solve the problem in curvilinear 
coordinates of the general form from the beginning. From a practical point 
of view the second approach is more convenient although it has been proved 
that this two approaches are mathematically equivalent. It should also be 
noted that its appear that the curvature of the boundary surface plays an 
important role in the boundary conditions obtained in both approaches. 

Solving the problem of clastic boundary conditions derivation, we get 
the opportunity to solve, in principle, any boundary problems needed to 
describe flexoelectricity in finite-size bodies. However, the corresponding 
boundary problems are extremely complex and difficult to solve. Such a 
boundary problem has been solved for a homogeneously polarized ball but 
for more complex geometries to make a similar is unrealistic. Even for the 
ball the solution is extremely cumbersome. Therefore the development of 
approximate methods for solving such boundary problems is desired. This is 
the second theoretical problem. 

The problem of developing a method for the approximate solution of the 
corresponding boundary problems were also solved and the solution is de¬ 
scribed above. It turns out that in the framework of such problems elastic 
displacements can be approximated as the sum of two parts. The first part 
was called non-classical, it is concentrated near the surface of the body and 
decays exponentially inside the body. To find this part it is necessary to 
solve fairly simple one-dimensional equations. The second, classical part is 
determined by the equations of the classical theory of elasticity. It turns 
out that the boundary conditions for this equations have a standard classical 
form of the boundary conditions for the body under the external forces on 
the surface. Thus, the finding of the classical part is reduced to the standard 
classical problem which does not require a separate discussion. It should 
be noted that these forces on the surface are formal. Physically there is 
no forces on surface, formal forces describe the interaction between classi¬ 
cal and non-classical part of the elastic displacement. So that these formal 
forces are expressed in terms of non-classical part of the displacement. The 
corresponding equations are given. 

It is important to note that for homogeneous polarization the formal 
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forces, describing the influence of the non-classical part of displacement to 
the classical one, appear only on the curved parts of the surface of the body. 
Thus, the above statement that curvature of the surface plays an impor¬ 
tant role gets a clear physical meaning. Note that the body can be, say, a 
polyhedron. In this case the formal forces appear only on the edges which 
can be treated as a limiting case of a curved surface. While calculating one 
should slightly smooth these edges and tend to zero the radius of smoothing 
at the end of the calculations. The similar approach is applied above for the 
particular case of the rod and thin plate which also has the sharp edges. 

Thus, in framework of the theory of continuum the description flexoelec- 
tricity in fini te-size bodies actually requires special theory. This theory is 
described above in detail. The application of this theory is illustrated in par¬ 
ticular problems of calculation of the bending of homogeneously polarized 
rod and a thin plate. 
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